arXiv: 1501.04991 vl [hep-ex] 20 Jan 2015 


Recent Results From The Daya Bay Experiment 

Chao Zhang 

(on behalf of the Daya Bay collaboration) 

Brookhaven National Laboratory, Upton, New York, USA 

Abstract. The Daya Bay reactor neutrino experiment has observed the disappearance of electron antineutrinos from nuclear 
reactors at —kilometer baselines. The relative measurement of the v e rate and spectrum between near and far detectors allows 
for a precision measurement of the oscillation parameters sin 2 20[3 and Am 2 c |. Two new antineutrino detectors (ADs) were 
installed in summer 2012, bringing the experiment to the final 8-AD configuration. With 621 days of data, Daya Bay has 
measured sin 2 20j3 = 0.084±0.005 and |Am 2 e | = 2.44+jJ'jj x 10 -3 eV 2 . This is the most precise measurement of sin 2 20!3 
to date and the most precise measurement of |A/M 2 f | in this channel. Several other analyses are also performed, including an 
independent measurement of sin 2 20j3 using v e samples tagged by neutron capture on hydrogen, a search for light sterile 
neutrinos, and a measurement of the absolute reactor antineutrino flux. 
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INTRODUCTION 

The discovery of neutrino oscillation over the past several decades has established that neutrinos have mass and their 
flavors ( V e , vy, V T ) mix. It represents one of the very few instances that indicate that the Standard Model is incomplete. 
The neutrino mixing can be characterized by three mixing angles (012,023,013). a possible CP violating phase S C p, 
and two mass-squared differences Atnji and Am^i. Determination of all their values with increasing accuracy was and 
continues to be the main goal of the experiments. 

The value of 0i3 has been a longstanding puzzle. Prior to 2012, only an upper limit of sin 2 2013 < 0.15 at 90% C.L. 
was obtained by the CHOOZ [1] and the PALO VERDE [2] experiment. The cleanest way to measure 0i3 is through 
kilometer-baseline reactor neutrino oscillation experiments. The reactor v e oscillation at ~km is dominated by the 0n 
terms, with the survival probability given by 

P ee = 1 — sin 2 2013(cos 2 012 sin 2 A31 +sin 2 0i2 sin 2 A32) — cos 4 0i3 sin 2 20i2sin 2 A21 (1) 

where A, ; - = 1.267Am 2 - Jfy'eV ) . Unlike the accelerator neutrino experiments, the reactor measurements are independent 
of 023 and the CP-violating phase Sep, and only slightly dependent on the neutrino mass hierarchy and the matter 
effect. In 2012, all three second-generation reactor ©13 experiments. Double Chooz, Daya Bay and RENO, reported 
clear evidences of v e disappearance [3, 4, 5] at --kilometer baselines with only a few month’s running. In particular, 
Daya Bay excluded 0j3 = 0 by 5.2 standard deviation with 55 days of data and provided the most precise measurement 
of 0i3. In 2013, Daya Bay updated the results with 217 days of data, with a spectral analysis to measure the oscillation 
frequency, which led to the first direct measurement of the v e mass-squared difference Am 2 e | 1 [6], In this article, 
we report the new results from Daya Bay with 621 days of data and two newly installed antineutrino detectors. It 
represents four times more exposure than the previously reported results in Ref [6]. 


THE DAYA BAY REACTOR NEUTRINO EXPERIMENT 

The Daya Bay reactor neutrino experiment was designed to provide the most precise measurement of 0; 3 among the 
existing and near future experiments, with sensitivity to sin 2 2013 < 0.01 at 90% C.L. This is achieved through a careful 


1 The effective mass-squared difference Am),, is defined such that sin 2 A ee = cos 2 0| 7 sin 2 A 31 + sin 2 0b sin 2 A 32 . 



site choice and detector design. A detailed description of the Daya Bay experiment can be found in Ref [7, 8]. 

The Daya Bay experiment is located near the Daya Bay nuclear power plant (NPP) in southern China. The six 
reactors, two at Daya Bay and four at Ling Ao, provide a total of 17.6 GW thermal power, making it one of the 
largest in the world. Prior to Aug 2012, Daya Bay had six functionally identical antineutrino detectors (ADs), each 
with a 20-ton gadolinium-doped liquid scintillating (Gd-LS) target region. A near-far arrangement of the ADs in three 
experimental halls (EHs) is used. Two ADs were placed in EH1 at a distance of 364 m from the two Daya Bay reactor 
cores. One AD was place in EH2 at a distance of about 500 m from the four Ling Ao reactor cores. And three ADs 
were places in EH3 1912 m from the Daya Bay cores and 1540 m from the Ling Ao cores. The relative measurement 
of v e rate and spectrum between near and far ADs allows for the cancellation of the relatively large uncertainty (2-5%) 
in predicting the absolute reactor v e flux. The functionally identical design of the ADs further removes the correlated 
detector systematics. 

Each AD consists of a cylindrical, 5-m diameter stainless steel vessel (SSV) that houses two nested acrylic 
cylindrical vessels. A 3-m diameter inner acrylic vessel (IAV) holds 20 tons of Gd-LS as the detection target. It is 
surrounded by a region with 20 tons of liquid scintillator (LS) inside a 4-m diameter outer acrylic vessel (OAV). 
Between the SSV and OAV, 37 tons of mineral oil (MO) shields the LS and Gd-LS from external radioactivity. Each 
AD has 192 8-inch PMTs mounted on the side walls of the SSV. Optical reflectors are place at the top and bottom of the 
SSV, which effectively increases the photo-covera ge to 12% . On average, 160 photoelectrons are collected per MeV, 
which results in an energy resolution of ^8%/y / £(MeV). Each AD has three automated calibration units (ACUs) 
mounted on top of the SSV lid. Weekly calibration is performed through remote deployment of an LED diffuser ball, a 
68 Ge source, and a combined source of 241 Am 13 C and 60 Co into the Gd-LS and LS liquid volumes along three vertical 
axes. The position accuracy of the deployment is < 7 mm [9]. 

The ADs in each EH are submerged in a water pool, with at least 2.5 m of high-purity water in all directions to shield 
against ambient radiation. Each water pool is optically separated with Tyvek sheets into inner and outer water shields 
(IWS and OWS) and instrumented with PMTs to function as Cherenkov-radiation detectors. The detection efficiency 
for long-track muons is > 99.7% [10]. The water pool is covered with an array of resistive plate chambers (RPC) for 
independent muon tagging. 

Antineutrinos are detected via the inverse beta decay (IBD) interaction, v e + p —> e + + n. The coincidence of the 
prompt signal from the e + and the delayed signal from neutron capture on Gd provides a distinctive v e signature. It 
is crucial to characterize the detector response between ADs due to the relative measurement strategy of the Daya 
Bay experiment. Besides regular calibration with ACU sources ( 68 Ge, 60 Co, 241 Am 13 C), special calibration sources 
( 137 Cs, 54 Mn, 40 K, 241 Am 9 Be and Pu- 13 C) were used in summer 2012 with ADI and AD2 in EH1. Furthermore, a 
manual calibration system (MCS) was installed on ADI during the special calibration period. The MCS allows for a 
full-volume deployment of sources to study the position dependence of the energy reconstruction. In addition, gamma 
and alpha peaks which could be identified in data ( 40 K, 208 T1, n capture on H, C, and Fe, 212 Po, 214 Po, 215 Po) are 
included. The difference in reconstructed energy between ADs is less than 0.2% in the energy range of reactor v e . 

Interpretation of the observed prompt energy spectra requires precise characterization of the absolute energy scale 
for e + , e~, and y. There are two major sources of nonlinearity between the true energy of the particle and the 
reconstructed energy: scintillator nonlinearity and electronics nonlinearity. The scintillator nonlinearity is particle 
and energy dependent, and is related to intrinsic scintillator quenching and Cherenkov light emission. It is modeled 
based on Birks’ formula and constrained by a combined fit with mono-energetic gamma peaks and the 12 B [3-decay 
spectrum, where the 12 B is a spallation product after muon interaction with carbon. The electronics nonlinearity is 
introduced due to the interaction of the scintillation light time profile and the charge collection of the front-end 
electronics. It is modeled with Monte Carlo and with the measurements from special single-channel FADC runs. 
The combined energy nonlinearity is cross-validated with internal radioactive 214 Bi, 212 Bi and 212 T1 j3-decay spectra, 
Michel electron spectrum in the ADs, and standalone bench-top Compton scattering measurements. The absolute 
energy scale uncertainty, which is correlated among all ADs, is constrained to be less than 1% in the majority of the 
energy range of reactor v e . 

Daya Bay took data in the aforementioned 6-AD configuration from December 2011 to July 2012. All previous 
results [5, 6, 8] were based on the data in the 6-AD period. In summer 2012, two additional ADs were installed in EH2 
and EH3, which completed the final 8-AD configuration of the experiment. The data taking resumed after October 
2012. The new results are based on the complete data set of the 6-AD period with the addition of the 8-AD period 
from October 2012 to November 2013, a total of 621 days. 



THE NEUTRINO OSCILLATION ANALYSIS 


IBD candidates are selected with the following criteria. First, events caused by PMT light emission are efficiently 
rejected. Candidates are then selected by requiring a prompt-like signal (0.7-12 MeV) in coincidence with a delayed- 
like signal (6-12 MeV) separated by 1-200 /is. Candidates are removed if the delayed-like signal occurs (i) within 
0.6 ms from a water pool muon, (ii) within 1 ms from an AD muon (7/ u > 20 MeV), or (iii) within 1 s from an 
AD showering muon (£), > 2.5 GeV). Finally, a multiplicity cut is performed to select only isolated candidate pairs. 
The combined efficiency in selecting v e Gd-capture events is 80.6% ±2.1%. The relative uncertainty in the efficiency 
between ADs is 0.2%, which is the dominant factor in determining oscillation parameters. A total of 1.1 million (150k) 
IBD candidates are selected in the near (far) halls, representing the largest v e sample among all previous and current 
reactor neutrino experiments. The time-dependent signal rates are strongly correlated with reactor v e flux, which is 
predicted from the thermal power, fuel burn-up, and exchange and enrichment records provided by the NPP. 

Daya Bay is a low-background experiment. The largest background comes from accidental coincidence of singles 
events, which accounts for 2.3% (1.4%) of the candidates in the far (near) hall. The accidental background, however, 
can be statistically calculated with high precision. A second background comes from neutrons from the ~0.7 Hz Am- 
C calibration sources inside the ACUs on top of the ADs. They can occasionally mimic IBD events by inelastically 
scattering with nuclei in the shielding material and then capturing on Fe/Cr/Mn/Ni, producing two y-rays that both 
enter the scintillating region. This background, which used to be the second largest in Daya Bay, is largely reduced 
after summer 2012 when the Am-C sources were removed from the two off-central-axis ACUs for the ADs in the far 
hall. It accounts for 0.2% (0.03%) of the far (near) hall candidates. Cosmogenic background in Day Bay is low due 
to good overburden of each underground hall (EH1: 250 m.w.e; EH2: 265 m.w.e; EH3: 860 m.w.e). The background 
from beta delayed-neutron emitters 9 Li and 8 He is 0.4% (0.4%) in the far (near) hall ADs. The fast neutrons produced 
by untagged muons account for 0.1% (0.1%) of the far (near) hall candidates. Finally, the background due to ( a,n ) 
nuclear reaction is 0.1% (0.01%) in the far (near) hall ADs. 

The left plot of Figure 1 shows the ratio of the detected to expected v e signals at the 8 ADs located in the three 
EHs, as a function of effective baseline. Consistent signal rates are measured for ADs in the same EH. The signal rate 
at the far site shows a clear ~6% deficit with respect to the near sites. The right plot of Figure 1 shows the measured 
background-subtracted spectrum at the far site compared to the expected spectrum based on the near site data both 
without oscillation and with the best-fit oscillation included. In the bottom panel, the ratio of the far site spectrum 
to the weighted near site spectrum is shown. The rate deficit and spectrum distortion are highly consistent with the 
oscillation interpretation. 

The oscillation parameters are extracted from a purely relative comparison between data at the near and far 
halls. The observed prompt energy spectra of the near halls are extrapolated to the far hall and compared with 
observation. This process is done independently for each prompt energy bin. A covariance matrix, generated from 
a large Monte Carlo data set incorporating both statistical and systematic variations, is used to account for all 
uncertainties. An alternative analysis using nuisance-parameter-based '/} method yields consistent results. Figure 2 
shows the allowed regions for sin 2 20i 3 and |A m 2 ee | at the 68.3%, 95.5%, and 99.7% confidence level. The best fit values 
are sin 2 20 13 = 0.084±0.005 and \Am 2 ee \ = 2.44+°;}° x 1(T 3 eV 2 with J 2 /NDF = 134.7/146. Although last known, 
the precision in the 0j 3 measurement (6%) is now the best among all three mixing angles. The |Am 2 e | measurement 
is highly consistent with and of comparable precision to the muon neutrino disappearance experiments. Under the 
assumption of normal (inverted) neutrino mass hierarchy, this result is equivalent to |Am 2 2 | = 2.39//]' }° x 10~ 3 eV 2 
(|A»4| = 2.49+U;}? x 1(T 3 eV 2 ). 

An independent measurement of sin 2 20 33 is performed using the V e IBD samples tagged by neutron capture on 
hydrogen (nH). Inside the Gd-LS region, ~15% of neutrons capture on H instead of Gd. In addition, the full 20-ton 
LS region is used as detection target. Due to the high accidental background caused by the long neutron capture time 
(~200 /is) and low y-ray energy (2.2 MeV) after n capture, the prompt candidate is required to have E > 1.5 MeV 
and the distance between the prompt and delayed candidates is required to be < 0.5 m. The analysis uses the complete 
217-day data set of the 6-AD period. A clear rate deficit is observed at the far hall. The rate-only analysis yields 
sin 2 20| 3 = 0.083 ±0.018, which is highly consistent with the results from nGd analysis. Details of the nH analysis 
can be found in Ref. [11], 

If a fourth (sterile) neutrino exists, it could cause additional spectral distortion betweens the ADs at different base¬ 
lines. An analysis to search for such sterile neutrinos is performed using a minimal extension of the standard model: 
the 3(active) + l(sterile) neutrino mixing model. Since Daya Bay has multiple baselines, the search depends only 
on the relative spectral difference between three sites and is largely independent of the reactor related uncertainties. 





FIGURE 1 . (left) Ratio of the detected to expected v e signals at the 8 antineutrino detectors (ADs) located in three experimental 
halls as a function of effective baseline. The expected signal is corrected with the best-fit reactor v e flux normalization. The 
oscillation survival probability at the best-fit value is given by the red curve, (right) The top panel shows the measured background- 
subtracted spectrum at the far site compared to the expected spectrum based on the near site data both without oscillation and with 
the best-fit oscillation included. The bottom panel shows the ratio of the far site spectrum to the weighted near site spectrum. The 
red curve shows the expectation at the best-fit oscillation values from the rate and spectral analysis. 



FIGURE 2. Allowed regions for sin 2 20)3 and |Am 2 e | at the 68.3%, 95.5%, and 99.7% confidence level, obtained from com¬ 
parison of the rates and prompt energy spectra measured by the near and far site antineutrino detectors. The best estimate of the 
oscillation parameters is given by the black dot. The adjoining panels show the dependence of A% 2 on each of the parameters. 


The analysis uses the complete 217-day data set of the 6-AD period. The relative spectral distortion due to the dis¬ 
appearance of v e is found to be consistent with that of the three-flavor oscillation model. The exclusion contours for 
sin 2 2014 and |Am 2 J are determined using both the Feldman-Cousins method and the CLs method [12]. The derived 
limits cover the l(r eV 2 < |Aw 2 j| < 0.3 eV 2 region, which was previously largely unexplored. Details of the sterile 
neutrino analysis can be found in Ref. [13]. 





































MEASUREMENT OF THE ABSOLUTE REACTOR ANTINEUTRINO FLUX 


The large reactor v e sample collected at Daya Bay allows for a precise measurement of the absolute reactor antineutrino 
flux. The analysis uses the complete 217-day data set of the 6-AD period. A total of 300k (40k) candidates are detected 
at the near (far) halls. Figure 3 shows the measured reactor v e event rate at each AD after correcting for the v e survival 
probability, re-expressed as Yo (cm 2 GW -1 day -1 ) and ay (cm 2 fission -1 ). The measurement among ADs is consistent 
within statistical fluctuations after correcting for the difference in the effective fission fractions. The uncertainty (2.3%) 
of the measurement is dominated by the uncertainty in detection efficiency (2.1%), which is correlated among all ADs. 
The measurement yields an average Yo = 1.553 x 10 -18 cm 2 GW -1 day -1 and oy = 5.934 x 10 -43 cm 2 fission -1 , 
with the average fission fractions 235 U : 238 U : 239 Pu : 241 Pu = 0.586 : 0.076 :0.288 : 0.050. 



FIGURE 3. The measured reactor v e event rate at each AD after correcting for the v e survival probability, re-expressed as Yq 
( cm 2 GW -1 day -1 ) and ay (cm 2 fission -1 ). The solid and open circles show the data without and with correction for the difference 
in the effective fission fractions observed by each AD. The uncertainty of the measurement is shown as the gray band. Three 
theoretical model predictions are shown as a reference. 

Three theoretical model predictions are shown in Figure 3 as a reference. The Huber [14] and ILL [15, 16] models 
predict the v e spectra for 235 U, 239 Pu and 241 Pu, while the Mueller [17] and Vogel [18] models predict for 238 U. The 
uncertainty in the model predictions is ~2.7%. The ratio (R) of the Daya Bay measurement to the Huber+Muller 
model prediction is R = 0.947 ±0.022, while R = 0.992 ±0.023 when compared to the ILL+Vogel model prediction. 



FIGURE 4. The reactor v e interaction rate of the 21 previous short-baseline experiments [19, 20] as a function of the distance 
from the reactor, normalized to the Huber+Mueller model prediction. [14, 17]. Experiments at the same baseline are combined 
together for clarity. The Daya Bay experiment is placed at the effective baseline of 573 m. The rate is corrected by the v e survival 
probability at the distance of each experiment, assuming standard three-neutrino oscillation. The horizontal bar (blue) represents 
the global average and its la uncertainty. The 2.7% reactor flux uncertainty is shown as a band around unity. 

A comparison between the Daya Bay result and the 21 past reactor neutrino flux measurements is performed 
following the methods in Ref. [19, 20]. The result of the comparison is shown in Figure 4. The Huber+Mueller model 
is used as the common reference model for all experiments. A neutron lifetime of 880.1 s [21] is used. The v e survival 




















probability is calculated with sin 2 26\ ] = 0.089 ± 0.009 determined from the rate-only analysis using the complete 6- 
AD data set [6]. The global average of the 21 past measurements with respect to the Huber+Mueller model prediction 
is determined to be R = 0.943 ±0.008 (experimental uncertainty), which is consistent with Rdyb = 0.947 ± 0.022 from 
the Daya Bay measurement. 


SUMMARY 

The Daya Bay experiment uses the relative measurement of the v e rate and spectrum between near and far detectors 
to precisely measure the oscillation parameters sin 2 20n and |Am 2 e |. Two new antineutrino detectors (ADs) were 
installed in summer 2012, bringing the experiment to the final 8-AD configuration. With 621 days of datas, Daya Bay 
has measured sin 2 20i3 = 0.084±0.005 and \/Xm 2 ee \ = 2.44+q}j x 10 3 eV 2 . This is the most precise measurement 
of sin 2 2013 to date. The precision measurement of 0i3 opens the door for future experiments to study neutrino mass 
hierarchy and leptonic CP violation. The \Am 2 e \ measurement is highly consistent with and of comparable precision 
to the muon neutrino disappearance experiments. By the end of 2017, Daya Bay expects to measure both sin 2 2013 and 
\Am 2 e | to precisions below 3%. 

Several other analyses are also performed. An independent analysis using the v e IBD samples tagged by neutron 
capture on hydrogen has measured sin 2 20i3 = 0.083 ±0.018. A sterile neutrino search has set stringent limits in the 
10 3 eV 2 < (Am 2 ! | < 0.3 eV 2 region. Finally, the absolute reactor antineutrino flux measurement has yielded consistent 
results with previous short-baseline reactor neutrino experiments. 
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